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1 Introduction 

Membrane-spanning peptides such as the ion channel gramicidin dimer, cause 
membrane deformation which contribute significantly to both the energetics of 
the insertion and the membrane-mediated interaction between the channels. 
Recently, interest in this field shifted from study of individual channels to study 
of collective effects in channel kinetics. Here we discuss an efficient way to 
describe the membrane-mediated interaction between the insertions in terms of 
coupled harmonic oscillators and corresponding " spring constants" . 

2 Variational principle, boundary conditions and 
Euler-Lagrange equation. 

We consider an elastic system extended in x, y— plane, where the deforma- 
tion can be described by a 2— dimensional field of "vertical" displacement u(r) 
(r = (x. y) is the radius vector in the mid-plane of the system). The examples 
include smectic and similar models for lipid bilayers , and " floating plate" model 
of classical elastic theory (see for review). 

The elastic boundary problem can be formulated as a variational (minimum) 
principle for the energy functional 

F^[u] = J gW{u,Vu,Au,...)df (1) 

where , the surface density of the elastic energy, is a quadratic function 
of the surface displacement u and its derivatives. We will consider a membrane 
with N cylindrical insertions, assuming that on the contour j i oli — th insertion 
both u(r) and V«(r) are fixed functions of r 7 . , the position vector for the points 
belonging to 7 i . It leads to the boundary conditions 

u ( r )k, = Ui(r 7 .) (2) 
Vu(r)\» =SiK) (3) 
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" n" designates the direction normal to r y i at the point r 7 . .Note that the ver- 
tical displacement Ui in the immediate contact of a membrane with an inserted 
peptide is typically described by the " hydrophobic matching condition" ||, [TJ 
leading to a particular case of Eq. |^ with Ui(r 7 .) = uq = const. 

Additional conditions on the external membrane boundary (designated as 
7oo) are: 

^(r)| 7oo =0 (4) 

Vu(r)\^=0 (5) 

In many cases, including the biharmonic problem with smooth and continu- 
ous boundaries that we consider ([Q, Ch. 6), the variational principle SF^ — 
(the minimum condition for the energy functional) leads to the Euler-Lagrange 
equation which we present for now as 

L(u) = (6) 

where L is a linear differential operator Q The elastic energy E = min III 
is the value of F^ 2 [u] calculated with the solutions of Eqs. ^-|| in place of u. We 
will show now that in some important cases E can be explicitly presented as a 
quadratic form of the boundary parameters, such as Ui and s^. 



3 Effective spring constants 



3.1 Boundary displacements and contact slopes fixed to 
constants 

Suppose now that the boundary displacements and the contact slopes are fixed 
at the i — th insertion to the constants 

Ui and s.j. Some preliminary results for this case were reported in Eqs. 
||, H can be written as 

u(r)| 74 =«i (7) 

Vu(r)|» = s t (8) 

We now introduce the " superfinite" elements, ^"(r) and <fii(r), solutions of 
Eq. H satisfying boundary conditions Eq. ||H and following conditions at the 
internal boundaries: 

ff(r)| 7fc = 8*, V#(r)|» fc = (9) 

1 Which means that L ■ (ci u\ + C2 U2) = ci L ■ u\ + C2 L ■ us 

(ci and C2 are the arbitrary constants). The linearity does not impose, however, any restric- 
tions on the order of the differential equation (which for most of the applications considered 
is biquadratic). 
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<^(r)| 7fc =0, V^(r)|^=<y« (10) 

where 5^ is the Kronecker symbol. The solution of Eqs. || -|], [?], || can be 
written as a linear combination of the superfinite elements: 

N 

«(r)=£(«i#(r)+*#(*)) (11) 

i=l 

Substituting this result into Eq. [l] allows presenting £ as a quadratic form 
of the boundary parameters: 

N N 

t=i j=i 

where c™' 3 are the effective spring constants describing interaction between 

the insertions i and j (c°f corresponds to the "self energy" of the i— th inser- 
tions). Such a "linear spring model" was first introduced for a single insertion 
in H and later generalized in Q . Eq. |l2| implies that the additional summation 
is performed over the repeated indexes, a, [3 (= u, s). 

The effective spring constants satisfy the symmetry relation 

c if = c f^ and ' consequently, cf? = c^ a . (13) 

As an illustration, consider the expression for g( 2 > which is typical in study 
of membranes: 

= B (Au) 2 + A u 2 (14) 

where A and B are proportional respectively to the membrane stretching 
and bending elastic constants and can be dependent on r, but not on u; A = 
d 2 jdx 2 + d 2 /dy 2 is the Laplace operator. 

Combining Eqs. [llj ^ and [l|, we find the spring constants: 

cf = (2 - Sy) J (BAtfA^ + A^)df (15) 
3.2 Azimuthal variation of the contact slope 
3.2.1 General formulas 

A possibility that the contact slope can become anisotropic at the contours of 
two interacting insertions was studied in |^|. The slope was presented as 

Si(r 7 .) = en + &i/i(r 7 .) (16) 
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where cij and fcj do not depend on r, and /i(r) are the fixed functions. 
In practice, functions /j(r) approximate the azimuthal variation of the contact 
slope. If energy is minimized over the slope parameters a, and bi , then different 
choices of /, (r) result in different families of trial functions s, (r 7 . ) . 

The surface displacement u(r) can be expressed now through the boundary 
parameters aj, 6, and Uj. 



i=l 



u(r) =y)[« i ^(r)+a i 0?(r)+6 i ^(r)] (17) 



where </>"(r) and <£"(r) satisfy respectively the boundary conditions || and 
|l0| , while </>j(r) satisfies the following conditions: 

0j(r)| 7h =O, V$(r)|^ =Wi(r 7l ) (18) 



The energy is still described by Eq. 12 with the spring constants defined by 



Eq. |l2j; the additional summation is now performed over the repeated indexes 
a,(3 — (u, a, b) . 



3.2.2 Applications for two insertions 



We consider now two identical insertions. Due to possible fluctuations, param- 
eters Ui, a,i and bi, and functions /i(r) for two insertions can still be different. 
The energy of two insertions is 



E = c53*(u?+u!) + c??(a5 + a!) + cg(&? + 6§) + 2c^ {u iai + u 2 a 2 ) + 

2c$(uxbi + u 2 b 2 ) + 2cf 1 {a 1 b 1 + a 2 b 2 ) + c^mu 2 + c?f cna 2 + cf 2 bib 2 + 
c"2 (uia 2 + w 2 ai) + c^( Ul b 2 + u 2 bx) + c°^{a x b 2 + a 2 h) (19) 



where we used the symmetry conditions c"f = c 22 , = = . 

We now consider equilibrium setting. This case was discussed in ||. The 
allowed functions /i(r) should satisfy the condition /i(r) = f 2 (—r), where r =0 
designates the midpoint between the insertions, with boundary parameters for 
both channels identical: u\ — u = u, a\ = a 2 = a, b\ = b 2 = b. Then, the 
energy can be presented as 

E = C uu u 2 + C aa a 2 + C bb b 2 + C ua ua + C ub ub + C ab ab (20) 

where 

Cafi = (2-6 aP )(2cg + <%£) (21) 
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Thus the total number of the effective spring constants is reduced to six. We 
intent to study the optimized ("relaxed", "equilibrium") slope, so that 

dE dE 

These conditions lead to 

a = — u, b = — u, A = C ab - AC aa Cbb (23) 

As a result, the total energy minimized over a and b can be written as 

E min = Ku 2 (24) 

_ Cbb(Cua — 4C aa C tlM ) + C a b(C a bC uu — C ua C u b) + C aa C 2 b 

We can now see that six effective spring constants C a p define the equilibrium 
slope parameters a and b and the equilibrium energy E m i n . 

The effective spring constants can be found from Eq. [IB] . Sometimes, how- 
ever, it is more practical to use the energy values E[{u, a, b}] defined for different 
sets {u,a,b} of the boundary parameters. 

For every distance d between the insertions, the elastic energy E[{u, a, b}] can 
be calculated numerically [jp. Then, all six spring constants can be found. The 
following equations illustrate this approach: 

_ £[{u,0,0}] _ £[{0,a,0}] r _ £[{Q,o, ft}] 
u z a z b z 

r _ E{{u,a,0}}~C uu u 2 ~C aa a 2 ^ _ E[{u, 0, b}} - C uu u 2 - C bb b 2 



C a 



h 



1 - LLU i 

ua uo 
E[{0,a,b}}-C aa a 2 -C bb b 2 



With these constants, the equilibrium (or minimized over the slope parameters) 
interaction energy profile E m i n (d) can be determined from Eq. EG. 



4 Final remarks 

It was shown that interaction energy between the insertions can be described 
in terms of effective spring constants accounting for the coupling between vari- 
ous degrees of freedom introduced through the boundary conditions. After the 
spring constants are defined, the equilibrium slope (which can in general become 
anisotropic) and corresponding interaction energy can be defined analytically. 
This approach is much more efficient than the direct energy minimization used 
in H. Some applications of this approach will be considered elsewhere. 
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